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Abstract 

In this paper we present the novel method for the generation of periodic em- 
bedded surfaces of nonpositive Gaussian curvature. The structures are related to 
the local minima of the scalar order parameter Landau-Ginzburg hamiltonan for 
microemulsions. The method is used to generate six unknown surfaces of Ia3d sym- 
metry (gyroid) of genus 21, 53, 69, 109, 141 and 157 per unit cell. All of them 
but that of genus 21 are most likely the minimal surfaces. Schoen-Luzzati gyroid 
minimal surface of genus 5 (per unit cell) is also obtained. 

PACS numbers: 61.20.-p, 64.75. +g, 68.10. -m, 02.40. -k 
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The surfactant molecules, which are the main ingredient of soaps and detergents 
have the ability to solubilize oil in water, two liquids which in the binary mixture 
at normal conditions are immiscible. This ability stems from their chemical struc- 
ture; a surfactant molecule has the polar and nonpolar segments at two ends and 
thus is simultaneously hydrophobic and hydrophilic. The term amphiphilic (from 
Greek word: loving both) molecule is used since one end (polar) of the molecule 
is well solubilized in water while the other (nonpolar) in oil. Hence the molecule 
preferably stays at the oil-water interface, forming a monolayer. At high concen- 
tration of surfactant the physical interface made of these molecules orders, forming 
periodic structures of various symmetries. Similar behavior is observed in systems 
of biological molecules (lipids) which in water solutions self assembly into bilayers. 
In 1967/68 Luzzati et al^~"^ observed that the type of ordering in the lecithin-water 
and lipid-water systems cannot be described by the arrangement of simple sur- 
factant aggregates such as cylinders, planes or spheres. They observed the cubic 
bicontinuous phase of the Ia3d symmetry (gyroid) where the lipid bilayers formed 
a highly curved smooth (embedded) surface of the same symmetry. Such surface 
divides the volume into two disjoint subvolumes. The NMR, SAXS (small angle 
x-ray scattering) and surfactant concentraction measurements indicate that these 
surfaces closely resemble triply periodic minimal surfaces^"^ i.e. surfaces charac- 
terized by zero mean curvature at every point. The latter belongs to the broader 
class of periodic surfaces of nonpositive Gaussian curvature. Since the discovery of 
minimal periodic surfaces in 1865 by Schwarz only one periodic embedded gyroid 
surface of cubic symmetry and genus 5 (per unit cell) has been discovered and fully 
characterized^ . Apart from that six more of Pn3m and Im3m symmetry are known, 
all of low genus. Here we present the general method which can he used to generate 
periodic surfaces of nonpositive Gaussian curvature. We prove the efficiency of this 
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method by generating six new gyroid structures of genus 21, 53, 69, 109, 14 1 and 
157. All these surfaces but that of genus 21 are minimal; we have also generated 
the Schoen-Luzzati gyroid minimal surface of genus 5 per unit cell. 

We point out that the unique characterization of the periodic surfaces using 
standard x-ray scattering methods is very difficult if the motif in the unit cell is 
not precisely known^"^^'^^. Our theoretical method combined with experimental 
technique will be extremely helpful in future discoveries of new surfaces in biological 
systems. 

The surfaces of surfactant systems find application in the production of the 
mesoporous silicate systems, where in the synthesis process the ordered surfactant 
surface is used as a template for the three dimensional polimerization of silicate-^^'-^^. 
One obtains an ordered silicate pore system with the symmetry and geometry of the 
surfactant template. We believe that our method can also be used for the design of 
new mesoporous structures. 

Our method is based on the Landau-Ginzburg model which has been proposed 
by Teubner and Strey-^^ and Gompper and Schick-'^^'-'^^ on the basis of neutron 
scattering experiments performed on microemulsion (homogeneous ternary mix- 
ture of oil, water and surfactant) and later experiments and theory of their wetting 
properties-*^^'-*^^. The Landau-Ginzburg free energy functional has the following form: 



where ^(</)) = g2(l>^ — go- Here </>, the order parameter, has the interpretation of 
the normalized difference between oil and water concentrations; (72, go are positive 
constants and /o can be of either sign. The last term in Eq.(l) is the bulk free 
energy and describes the relative stability of the pure water phase (0 = —1), pure 
oil phase {(p — 1) and microemulsion [(p — 0). The stability of bulk microemulsion 
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phase depends on /q: for /o > microemulsion is a metastable bulk phase whereas 
pure water phase or pure oil phase are stable; for /o < microemulsion is stable. 
We note that in general g{(f)) can be a polynomial in 0^. 

For go > 2 the system can undergo a transition to periodically ordered phases 
where water rich domains and oil rich domains order. The interface between the 
domains corresponds to ^(r) = 0. The different structures (stable or metastable) 
correspond to the minima of the functional (1). The following simple argument 
shows that among the surfaces, inside these structures, we might expect minimal 
surfaces. The mean curvature of the surface at point r is given by the divergence 
of the vector normal to the surface at this point ^O-^^: 

\mj 2|v0i ^ ^ ^ 

Here Vn denotes the derivative along the normal to the surface. It follows from the 
second term of Eq(l) that F[(p] is minimized when |V^| has the maximal value for 
0(r) = since at that point g{4>) has the lowest value. For the maximum of | V0| its 
normal derivative vanishes and consequently the second term in Eq(2) does so. We 
also know that in the case of 4>, —(f) symmetry, H averaged over the whole surface 
should be zero. It means that either A0 is exactly zero at the surface or it changes 
sign. From the first term of Eq(l) it follows that the former can be favored and 
consequently H — Q aX every point at the surface. Hence we can expect that some 
of the surfaces are minimal. This argument does not take into account the global 
distribution of the field </>, nonetheless it provides a useful hint for our studies. 

In order to find the minima of the functional we have discretized Eq(l) on 
the cubic lattice. Thus the functional F[0(r)] becomes a function F{{(f)i^j^j^}) of 
variables, where Nh is the linear dimension of the cubic lattice and h is the 
distance between the lattice points. Each variable (/>i,j,fe represents the value of the 
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field 0(r) at the lattice site {i,j,k), and the indices i,j,k change from 1 to N. In 
our calculations we use N=17,33,65 and 129; final results are shown for N=129. 
Please note that N=129 results in over 2 milion points per unit cell. The first 
and second derivatives in the gradient and laplasian terms of the functional (1) 
were calculated on the lattice according to the three point formula for the first 
derivatives and five point formula for the second derivatives We impose on the 
field (f)ijf~ the periodic boundary conditions and the symmetry of the structure, we 
are looking for, by building up the field inside a unit cubic cell from a smaller 
polyhedron, replicating it by refiections and rotations combined with translations, 
since the gyroid symmetry involves glide planes. Such procedure enables substantial 
reduction of independent variables. 

The initial configuration needed for the minimization is set up by building the 
field 0(r) first on a small lattice A?" = 3 or 5. It is done by analogy to the structure 
of a two component (A,B) molecular crystal. The value of the field 4'i,j,k at a lattice 
site (i, j, k) is set to 1 if in the molecular crystal an atom A is in this place. It is set 
to —1, if there is an atom B, and to 0, if there is an empty place at (z, j, k). Next the 
small lattice can be enlarged to desired size by changing the number of points from 
N to 2A'" — 1 and finding the values of 0i,j,fc in new lattice sites by interpolation. 

We have used the conjugate gradient method^^ to find a minimum of the func- 
tion F{{(f)i,j^k})- It is highly unlikely, because of numerical accuracy, that a value 
of the field (f>i,j,k at a latice site (i, j, k) is exactly zero. Therefore the points of the 
surface have to be localized by linear interpolation between the neighbour sites of 
the lattice. This approximation is legible because the field 0(r) is very smooth. The 
points on the surface are used for the triangulation. From the triangles covering 
the surface we get the surface area and the Euler characteristic, %• The latter is 
given by the Euler relation: x = F — i? + where F is the number of faces, E is 
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the number of edges and V is the number of vertices of the triangles covering the 
surface. The edges and vertices has to be taken with weight 1, 1/2, or 1/4 if they 
appear inside, at the face, or at the edge of the unit ceU, respectively. 

We have performed the detailed study of the phase diagram, checking the 
Landau free energy of almost 30 different structures of various symmetries and 
found that the only stable ordered structure is the lamellar phase. The phase 
boundaries for the lamellar phase are given by Gompper and Zschocke^^. The 
gyroid phase of genus 5 has the second lowest energy (after lamellar phase) among 
all the structures (Table 1). We have noted that this phase has larger area per 
unit volume than the lamellar phase, thus in the case of very sharp interfaces it 
should have smaller energy than the lamellar phase, since from the second term 
of Eq(l) it follows that at 0(r) = the gradient term gives large and negative 
contribution to the energy. This negative contribution is not cancelled by the bulk 
or laplacian terms. Unfortunately the size of the interface scales with the size of 
the unit cell. We have performed the same calculations for the new function gi^cf)) 
given by : g{(j)) = ^20^ — ^o- Increasing the power of by a factor of two indeed 
sharpens the interface between oil and water but at the same time reduces the size 
of the unit cell. The net result is the larger relative difference in energies between 
the gyroid and lamellar phase. We note that in the case of multiparameter Landau 
models introduced in recent years^^'^^ we may expect the stabilization of the various 
phases which here are only metastable. 

Among the local minima of the functional (Eq(l)) we have found four known 
minimal surfaces: P, D, I-WP and G^^. Here we present the detailed study of the 
surfaces of gyroid symmetry. In all these structures a surface is characterized by non- 
positive (zero or negative) Gaussian curvature. In Table 1 the main characteristics 
of the gyroid surfaces are given. Here we have used = 129 points per edge of the 
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unit cell. In order to estimate the errors we have compared the results obtained for 
N = 65 and N = 129; for the surface area the largest errors (gyroid 141 and 157) are 
smaller than 0.3%, for the energy the largest error is smaller than 1% for the gyroid 
5 structure and few percent for other structures. Of course these are the upper limits 
and most probably the errors are much smaller. In Figs. (1,2) two gyroid structures 
are shown together with the histograms of their mean curvature. Please note that 
for the minimal surface the mean curvature is peaked around in the histograms, 
but due to the numerical accuracy the peak has a finite width. We note that as 
the genus of the surface increases the surface area per unit volume (S/d^) and the 
energy per unit volume does not change very much(Table 1). Some authors^^ ruled 
out the possibility of existance of high genus surfaces in real systems, because of 
expected high curvature regions. We have checked the Gaussian curvature in high 
genus surfaces and found that it is not much different from the low genus surfaces. 
This is due to the sufficiently large size of the unit cell for the former structures. 
We also observe that the gyroid structures of high genus are most easily generated 
(from any initial configuration and sufficiently large unit cell) close to the stability 
region of microemulsion. 

Summarizing: We have used the Landau-Ginzburg hamiltonian for microemul- 
sions to generate periodic surfaces of nonpositive Gaussian curvature. Using this 
model we have studied 7 periodic gyroid surfaces: one of them is the Schoen-Luzzati 
minimal gyroid surface of genus 5. The remaining six structures are new. The model 
can be well applied by physicists working in soft condensed matter, mathematicians 
working in topology, biologists and crystallographers. We are positive that its rich- 
ness is far from being explored by our work. 

This work was supported by the KBN grants 2P03B 01810 and 30302007 and 
Foundation FWPN. 
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Table Caption 

Table 1. The gyroid surfaces of nonpositive Gaussian curvature. The symmetry of all 
the structures is Ia3d. In all cases the volume fraction is 0.5 by construction. 
Here g2 = 4-\/l + /o + + 0.01, = 3 and /o = 0. At this point the energy 
of the stable phase (lamellar phase) is -0.2077 and the size of the unit cell is 
d = 3.4. In column 2 the energy per unit volume is given. In the third column 
the dimensionless linear size of the unit cell, d is given. The surface area, 
S, (fourth column) is divided by y^/^ = d"^ i.e. is calculated per face of the 
unit cubic cell. The surface area per unit volume, S/d^, is almost constant for 
all the structures. The genus, g, (fifth column) has been calculated from the 
formula 1 — x/2, where x is the Euler characteristic per unit cell. We give, in 
the last column, the quantity 5 = \x\^^^d'^ / S. This quantity characterizes not 
only the ordered phase, but also the fluctuating microemulsion^^. We think 
that it can be used as a test for the structure of microemulsion. We find 
these structures practically for all values of the parameters where the lamellar 
phase is also stable^^, although most easily they are generated close to the 
microemulsion stability region. The genus, surface area per side of the unit 
cell S/d'^, symmetry, volume fraction and 6 do not depend on the parameters, 
go,fo and g2- Only the energy, size of the unit cell and surface area per unit 
volume are model dependent. 

Figure Captions 

Fig.l G141 gyroid structure (see Table 1). (a) One unit cell (b) The histogram of 
the mean curvature H. The size of the unit cell is given in Tablel. Herethe 
typical curvature, 1/R — ±\/—K, is ±0.1 {K is the Gaussian curvature). 

Fig.2 GM157 gyroid structure. Legend as in Fig.l. Here the typical curvature is 
±0.5. 
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